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Abstract 

Through laboratory measurements, we compare the rotation of spherical and ellipsoidal par- 
ticles in homogeneous, isotropic turbulence. We find that the particles' angular velocity statistics 
are well described by an Ornstein-Uhlenbeck (OU) process. This theoretical model predicts that 
the Lagrangian autocovariance of particles' angular velocity will decay exponentially. We mea- 
sure the autocovariance using stereoscopic particle image velocimetry (SPIV) applied to particles 
whose size is within the inertial subrange of the ambient turbulence. The SPIV resolves the 
motion of points interior to the particles, from which we calculate the solid body rotation of 
the particles. This provides us with the angular velocity time series for individual particles. 
Through ensemble statistics, we determine the autocovariance of angular velocity and confirm 
that it matches the form predicted by an OU process. We can further use the autocovariance 
curve to quantify the turbulent rotational diffusivity. 



1 Introduction 



In this paper, we compare rotational dynamics of spherical and non-spherical particles suspended 
in a turbulent flow. Rotation of particles is important in environmental and industrial pro- 
cesses, such as floc c ulatio n and papermaking, in which rotation can affect particle aggregation 



fcoch and ShaofehJ . H). In these and many other applications, particles are typically non- 
spherical. This leads to close coupling of the rotational and translational motion s and further 
complicates the mathematical treatment ( Mortensen et al. . 2008bl ; Clift et al. . 20051 ). 

For both spherical and non-spherical particles, rotation significantly effects th e wake structure 
behind a particle and t hus the forces coupling the particle and surrounding fluid (|Batchelorl . 119671 ; 



Giacobcllo et 



al 



20091). These forces, in tu r n, set the parti c le rot ation rate by causing angular 
acceleration ( Bagchi and Balachandar , 20021 ; Shin and Koch , 2005 ). This feedba ck can lead to 

' 20041 ): 



complex behaviors, e.g. broken symmetry in quiescent fluid (jJennv et all 120041 ) : rotation in a 
direc tion opposite to local vort i city d uring the transient approach to equilibrium in a steady shear 
flow (|Bagchi and Balachandarl . 120021): and spherical particles that do not rotate at half the local 
fluid vorticity in turbulent flow ([Mortensen et al.1 . 12001 l2008ah . 

Given the complexities of directly predicting particle rotation in turbulent flow, a stochastic 
approach is a useful alternative. In analogy to a translational random walk, we will model the time- 
dependent evolution of a particles' angular velocity as a random walk in rotation space. A stochastic 
process that closely approximates the translational motion of fluid parcels in turbulence is the 
Ornstein-Uhlenbeck (OU) process. The properties of this process are described in the next section 
but, importantly, it describes a diffusion process that is stationary, Markovian, and Gaussian. 
This paper hypothesizes that the OU process can be applied to the angular velocity of spherical 
and ellipsoidal particles, resulting in rotational diffusion that is analogous to the translational 
diffusion. We test this hypothesis by measuring the Lagrangian autocovariance and probability 
density functions (PDFs) of particles' angular velocity. Understanding rotational diffusion will 
allow prediction of the dynamics of particle orientation in turbulent flows 
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2 Theory 



The Langevin equation is a stochastic differential equation that describes the Ornstein-Uhlenbeck 
(OU) stochastic process. The OU process is completely characteri zed by a Gau ssian distribution 
and an autocovariance, both of which are statistically stationary ( Doob . 19421 ) . This stochastic 
process is Markovian, meaning that the probability distr ibution of the next state is determined 
only by the probability distribution of the current state ( NorrisL 19981 ). In the OU process, the 
autocovariance takes the specific form of a negative exponential, e.g. R(t) = a 2 e~ T l T , where a 2 
is the variance of the Gaussian distribution, r is the time-lag, and T is the integral timescale. 
The OU proces s has been s hown to accurately describe the Lagrangian velocity of fluid parcels in 
turbulent flow (|Popd . l200d ). One important implication of this is that OU process and Langevin 
equation can be used to accurately model the correlated random walk that describes the diffusion 
of a cloud of particles. The time-dependent diffusivity is given by T(t) = J Q * R(r)dr, and can be 
used to apply Fick's law for mass flux in a turbulent flow. 

We hypothesize that the OU process and Langevin equation are also a good model for the 
rotation of large particles suspended in a turbulent flow. Specifically, we hypothesize that the 
Lagrangian timeseries of angular velocities will follow an OU process. We can test this hypothesis 
by determining whether the Lagrangian angular velocity autocovariance decays exponentially as 



afe 



(1) 



defining Tn i as the integral timescale of Lagrangian angular velocity, and af as the angular velocity 
variance. We also check whether the PDF of particles' angular velocity is stationary and Gaussian. 
In §J]we see that observations match these predict i ons q uite well. As a result, we can find the time 



dependent rotational turbulent diffusivity (jTavlorl . 11921 



as 



TnM= J Rn i {r)dr = a.iT Qi [l-e 



(2) 



The rotational diffusity allows us to calculate the characteristic time for a complete rotation cycle. 



3 Methods 

We study particle motion in homogeneous, isotropic turbulence generated in a laboratory facility. 
This is a symmetrically stirred rectangular tank driven by a planar forcing element on each end. The 
symmetry contributes to large scale isotropy, as does the use of two ran domly actuated synthetic 
jet arrays (RASJAs) for the forcing elements ( Variano and Cowen . 20081 ). Each RASJA is a planar 



array of 64 synthetic jets. The jets fire in a stochastic pattern that maximizes the shear production 
of turbulence, similar to an active grid wind tunnel but without mean flow. The RASJAs allow the 
tank-scale secondary circulations to be essentially eliminated, which increases the residence time 
of particles in the image area helping us measure Lagrangian particle trajectories. Two screens 
are placed between the pumps and the center of the tank, which keep the particles contained in 
the central me asurement reg i on. T he flow properties of the RASJAs without screens are discussed 
extensively in iBellani et all (|2012h . The turbulent statistics for the RASJAs with screens and 



without particles can be found in table [TJ Although the screens decrease the turbulent kinetic 
energy, they do not affect the homogeneity or isotropy of the flow. 

We measure the turbulent statistics of the fluid-phase (without particles) using 2D particle 
image velocimetry (2DPIV). In the center region of the tank, the Reynolds number (defined with 
respect to the Taylor microscale) is ReA = 110. Single-point turbulent statistics are homogeneous 
and isotropic throughout a large center region. In previous experiments we found that this region is 
much larger than the turbulent integral length scale; therefore, the measurements are not affected 
by spatial heterogeneity. 
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Transverse velocity rms, 


Urms (m s 1 ) 


Si 0.012 




Vertical velocity rms, 


i — 1\ 

Vrma (m S 1 ) 


0.012 


[0.012 0.012] 


Longitudinal velocity rms, 


Wrms (m S i ) 


0.013 


[0.013 0.013] 


Turbulent kinetic energy, 


k (= ^(Ei«t ms )) (m 2 s 2 ) 


2.3 x 10~ 4 


[2.3 x 10~ 4 2.3 x 10~ 4 ] 


Taylor microscale, 


\ z (m) 


0.3 X 10 


tin iwin — 3 1 n — 3l 
[O.ZX1U 1UX1U J 


Integral length-scale, 


A z (m) 


57 X 10~ 3 


[57 x 10" 3 58 x 10" 3 ] 


Eddy turnover time, 


T (= A z /w rms ) (s) 


4.3 


[4.3 4.3] 


Kinematic viscosity, 


v (m 2 s _1 ) 


1.0 x 10~ 6 




Reynolds number (Taylor), 


Re A (= w rms \ z /v) 


110 


[81 140] 


Reynolds number (Integral), 


ReA (= Wrms^z/v) 


760 


[750 760] 


Turbulent dissapation rate, 


e (= l<iu[w 2 rms /\l}) (m 2 s- 3 ) 


3.8 x 10~ 5 


[1.9 x 10~ 5 5.8 x 10~ 5 ] 


Kolmogorov length-scale, 


q (= (^ 3 A) 1/4 ) M 


0.40 x 10" 3 


[0.40 x 10" 3 0.40 x 10~ 3 ] 


Kolmogorov time-scale, 


r„ (= (i// e )Va) (s) 


0.16 


[0.16 0.17] 



Table 1: Turbulent statistics of fluid-phase without particles in experimental facility with 95% 
confidence intervals. Both A 2 and A z are calculated from the longitudinal autocovariance of fluid- 
phase velocity. The introduction of particles significan tly reduces the tur bulent kinetic energy; 
spheres reduce k by 15% and ellipsoids reduce k by 3% (jBellani et al.l . 12012 ). 



d s (m) 


8 x 10" 3 


d e (m) 


16 x 1(T 3 


l e (m) 


8 x 10~ 3 


P P (kg m" 3 ) 


1007 



Table 2: Quantitative description of hydrogel particles. 

We place hydrogel spheres or ellipsoids (made of an agarose solution of 4 g/L) into the stirred 
tank. The spherical particles have diameter of d s = 8 mm and the ellipsoidal particles have polar 
and equatorial axes of l e = 16 mm and d e = 8 mm respectively. These dimensions are such that they 
match the Taylor length scale (from table [lj X z = 8.3 mm) and their relaxation time scale (3.64 s for 
spheres and 5.46 s for ellipsoids, calculated using the major axis) is close to the eddy turnover time 
of the ambient turbulence ( T = 4.3 s). The particle R eynolds numbers are much greater than 1 



and are 22 and 63 for spheres and ellipsoids respectively (jBellani and Variand . 120121 ) . The particles 
are near-neutrally buoyant (1007 kg m -3 ), thus they stay suspended during experiments. They 
are refractively matched to water and nearly transparent. These optical properties are essential, as 
they allow us to pass a laser sheet through the particles and measure the motion of interior points 
via stereoscopic particle image velocimetry (SPIV). Interior points are made visible by reflective 
tracers embedded in the hydrogel. 

Using a SPIV system (commercially provided by LaVision, Inc.), we find the flow field for the 
fluid-phase turbulence as well as the velocity vectors inside the particles. Our tempora l resol ution 



is 14.773 Hz, which is faster than the particle relaxation timescale (see Bellani et al. . 20121 ) and 



our spatial resolution is 1.38 mm, which is smaller than the particle size. We measure velocity 
vectors inside the solid particles by tracking clusters of the embedded tracers via standard SPIV. 
The resulting three-component velocity measurements along a 2D plane inside each particle allows 
us to calculate the particle angular velocity, fl = [Q, x , Q y , Cl z ]. We do so by using the solid body 
rotation equation: 

Vi=V M |f!xR t „, (3) 

where is the velocity at any interior point L, Vm is the velocity at a second interior point 
M, fi is the rotation about the center of mass, and Rlm is the distance vector between points 
L and M. Due to the fact that we measure three components of velocity at each point (e.g. 
Vl = [Vl x , V Ly , Vl z }) and all points are coplanar (i.e. ~R LM = [Rlmx, Rlm v ,0}), this system of 
equations is over-determined in the z-coordinate and under-determined in the x- and y-coordinates. 
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By including a third interior point (N), equation (|3|) can be rewritten twice, once in terms of L 
and N as well as once in terms of M and N. Using equations ([3]) and the two other permutations, 
all three components of rotation can be determined. This method gives a single value for ft x as 
well as Q y and four estimates of fl z . We average the four estimates of fl z , therefore our result 
for fl z contains less measurement noise than the other two components. Since there are multiple 
velocity points inside of the particle, we calculate ft once for each possible triplet, Vl, Vn and Vm- 
Compiling all of the measurements of ft for a single particle gives a distribution for each component. 
The median value fo r each component is selected as the best estimate of particle angular velocity 



(jBellani et all 12012 ). 



We measure ft for a particle as it travels along a trajectory through the turbulent flow. These 
trajectories range in length from 8 to 119 SPIV frames. Collecting ft for all frames gives us a time 
series of particle rotation. We then compute the autocovariance for each component of fli(t) from 
its definition: 

i?n.(r) = {^(t + r)O i (i)). (4) 

Here (•) signifies an ensemble average across all particle trajectories. 

We collect sufficient trajectories to calculate statistically converged data in the first and second 
moments of the distribution of That is, adding more data would decrease the uncertainty on 
the moments, but not change their values. Our data show convergence for the mean and variance 
of ft after 300 independent time series. In total, we collected 407 ellipsoids trajectories and 572 
sphere trajectories. With this data, we were able to also compute a fourth standardized moment 
(kurtosis) that conv e rged for fl z (the low-noise component) but not for the other two components. 



Snvder and Lumlevl (|197lh predict that 700 independent time series are needed to reduce the error 



below 10% in the measurement of the autocorrelation of translational velocity statistics. The 
value of 700 independent time series is slightly larger than value of 300, which we measure. This 
makes sense because we measure the autocovariance instead of the autocorrelation, and the former 
converges with fewer samples. 

The definition of ft used here is the angular velocity of a particle about its center of mass, 
with the three rotation components referenced to a set of coordinate axes that are fixed in the 
laboratory frame. It would also be possible, though more complicated experimentally, to decompose 
the rotation into the components relative to a set of coordinate axes that move with each particle. 
The two definitions will likely not give identical results, and both would be interesting to know. 
For example, when studying organisms, the rotation relative to the local axes would be useful 
in quantifying changes to motile diffusion and drag, and the rotation relative to the laboratory 
axes would be useful in determining the statistics of the organism's sensory search of 3D space. 
Herein, we focus only on the rotation relative to the lab coordinates, and evaluate this rotation 
while following each particle on its translational random walk. This is the exact meaning of our 
definition of "Lagrangian angular velocity," and we caution readers to be aware of the other possible 
definition. 

Because our data is transformed several times between the raw measurement (tracer locations 
in a plane) to the final measurement (particle angular velocity vector), it is nontrivial to propogate 
error to find the uncertainty in our measurements. Thus we use a Monte Carlo (MC) simulation 
to evaluate the uncert ainty of ft . We b egin by constructing synthetic angular velocity trajectories 
using an OU process (|Gillespiel . ll9Q6h . We then construct corresponding planar velocity vector 
fields of the type which we measure using SPIV. To these vector fields data, we add random noise 
that is normally distributed with zero mean. Finally, both the noisy and noise-free simulated vector 
fields are passed through the rotation algorithm to calculate ft. By comparing the results from the 
noise-free and noisy simulated measurements, we can assess the effects of measurement noise on 
the statistics of ft; the results are discussed in the next section. 
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Spheres 


95% CI 


Ellipsoids 


95% CI 




(Rad s" 1 ) 


0.43 


[0.39 0.47] 


0.57 


[0.54 0.61] 


a y 


(Rad s- 1 ) 


0.50 


[0.40 0.60] 


0.42 


[0.39 0.46] 


a z 


(Rad s" 1 ) 


0.47 


[0.44 0.49] 


0.39 


[0.38 0.41] 


Tq, 


. (8) 


0.43 


[0.29 0.57] 


0.37 


[0.30 0.43] 




, ( S ) 


0.48 


[0.15 0.81] 


0.37 


[0.28 0.45] 


To, 


, (s) 


0.39 


[0.33 0.45] 


0.52 


[0.45 0.59] 



Table 3: Parameters from the exponential fit: noise-free standard deviation (ctj) and Lagrangian 
integral timescale (ThJ. 

4 Results 

The experimental autocovariance curves for the spherical and ellipsoidal particles are shown in 
figures and [lb. Equation (pQ) is fit to the data with free parameters af and Th i5 reported in table 
[3j The characteristic timescale (TqJ of exponential decay found in these fits is the Lagrangian 
integral timescale of particle rotation, and af is the variance of particle angular velocity . 

Evident in figures [T] and is a considerable deviation between data and model at zero lag (r = 0) 
due to measurement noise. If this measurement noise is uncorrelated with itself over time, then the 
noise will be present only in the first point on each curve. Until lags greater than about a second 
and a half, this appears to be the case, because the data and exponential curve fit agree well. For 
larger lags, because there are many more short trajectories than long trajectories in our dataset, 
the fit and data are more discrepant (see figure [2]). The exponential fit allows us to estimate the 
noise- free variance at Rni( T = 0), which is equal to a?. 

We approximate the magnitude of measurement noise as the difference between the variance 
measured as the second moment of the angular velocity PDF (af) and the noise- free variance 
measured from the exponential fit to autocovariance {af). From these values we define a signal-to- 
noise ratio as 

SNR, = J (5) 

The SNR, measured variance, and noise-free variance are all reported in table 21 The SNR is an 
essential input into our Monte Carlo analysis; we set the simulated measurement noise in the Monte 
Carlo simulation so that its SNR matches that which we observe in the measurements. We then 
use the MC results to infer how measurement noise affects statistics other than the angular velocity 
variance. 

The rotation statistics measured from our data are given in table [3} along with corresponding 
MC results. Due to the isotropy of the turbulence in the tank (see table [1]), the rotation of the 
particles should also be isotropic. The fitted parameters on the exponential model (table[3|) indicate 
that this is true, to within measurement uncertainty. As expected, the effect of measurement noise 
is greater for rotation about the x- and y- axes than the z-axis. This is because in the algorithm to 
compute rotation, Vt z is the overdetermined component and therefore contains more data, increasing 
accuracy (see £J3J). This anisotropic noise effect can be seen clearly in the Monte Carlo results (table 
13]); when statistically isotropic rotation (seen in the column marked "Initial") is measured with our 
anisotropic imaging routine, the resulting PDF has anisotropic moments. 

Noise-free variance can be easily calculated as discussed above, but no such method is available 
for kurtosis. Thus to understand kurtosis we use the MC analysis as a guide. MC results show that 
our measurement method gives kurtosis values that are strongly upward-biased. Thus a distribution 
of rotation values that is truly Gaussian (K = 3) will appear super-Gaussian in our measurements. 
This bias is present in all coordinate directions, though it is stronger in x and y than in z. The 
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Figure 1: Autocovariance of rotation for (a) spheres, (b) ellipsoids. Top panel is £l x , middle panel 
is Vt y , and bottom panel is Q z . Points are experimental data. Solid line is exponential fit. Error 
bars are standard error. 
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Figure 2: Autocovariance of rotation about the z-axis for spheres (points and solid line) and 
ellipsoids (circles and dashed line). For both plots, error bars are standard error. Due to a large 
number of short trajectories, fit is closer for small lag time, r. 



biased measurements of kurtosis computed from our experiments are seen in table HJ only the z- 
component is statistically converged, thus we do not report confidence intervals for the others. The 
measured values of K z are close to the MC predictions for a biased measurement of a Gaussian 
distribution. In fact, they are slightly smaller than these predictions. We conclude from this that 
the true distribution of angular velocities in our experiment is near Gaussian, and may even be 
slightly sub-Gaussian (K < 3). 
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Stereoscopic Particle 


Image Vclocimctry (SPIV) 




Monte Carlo (MC) Simulation 






Spheres 


95% CI 


Ellipsoids 


95% CI 


Initial 


95% CI 


Pinal 


95% CI 






0.81 


[0.72 0.91| 


1.0 


[0.92 1.1J 


- 


- 


0.85 


[0.89 0.95| 


SNRy 




0.88 


[0.66 1.1] 


1.0 


[0.92 1.2] 


- 


- 


0.87 


[0.90 0.97] 


SNR 2 




1.4 


[1.2 1.5] 


1.1 


[1.1 1.2] 


- 




1.3 


[1.3 1.4] 


(Six) (Rad 




-0.012 


[-0.036 0.010] 


-0.024 


[-0.050 -0.00043] 


0.033 


[0.014 0.052] 


0.018 


[-0.014 0.047] 


{(l y ) (Rad 




-0.029 


[-0.054 0.0072] 


-0.052 


[-0.077 -0.032] 


0.033 


[0.014 0.052] 


0.032 


[0.0023 0.061] 


(Cl z ) (Rad 


B- 1 ) 


0.021 


[0.0042 0.036] 


0.011 


[-0.0041 0.025] 


0.033 


[0.014 0.052] 


0.031 


[0.0063 0.058] 


ff E (Rad s~ 


- 1 ) 


0.68 


[0.66 0.69] 


0.80 


[0.79 0.82] 


1.0 


[0.99 1.0] 


1.5 


[1.5 1.5] 


ay (Rad s~ 




0.76 


[0.73 0.79] 


0.59 


[0.57 0.61] 


1.0 


[0.99 1.0] 


1.5 


[1.4 1.5] 


a z (Rad s~ 




0.58 


[0.56 0.60] 


0.52 


[0.50 0.53] 


1.0 


[0.99 1.0] 


1.3 


[1.2 1.3] 


a x (Rad s~ 




0.43 


[0.39 0.47] 


0.57 


[0.54 0.61] 


1.0 


[0.99 1.0] 


0.97 


[0.94 1.0] 


CKy (Rad s" 




0.50 


[0.40 0.60] 


0.42 


[0.39 0.46] 


1.0 


[0.99 1.0] 


0.98 


[0.95 1.0] 


a z (Rad s 




0.47 


[0.44 0.49] 


0.39 


[0.38 0.41] 


1.0 


[0.99 1.0] 


1.01 


[0.99 1.0] 


Kx 




SJ2.7 




SJ2.7 




3.0 


[2.9 3.0] 


8.2 


[7.3 9.2] 


'<>, 




»5.1 




a;2.9 




3.0 


[2.9 3.0] 


7.7 


[6.8 8.7] 






5.7 


[4.9 6.7] 


5.3 


[4.8 5.9] 


3.0 


[2.9 3.0] 


6.6 


[6.0 7.2] 



Table 4: Moments of angular velocity PDF and Monte Carlo simulation PDF. Here Oi is the 
standard deviation of the PDF of f2j and K{ is the kurtosis of the PDF of f^. The SNR is calculated 
by measuring the noise magnitude as the difference between the directly measured standard devation 
(cij) and the noise-free standard devation (Qj) obtained from the exponential fit. 



Stereoscopic Particle Image Velocimetry (SPIV) Monte Carlo (MC) Simulation 





Spheres 


95% CI 


Ellipsoids 


95% CI 


Initial 


95% CI 


Final 


95% CI 


a x (Rad s~ x ) 


0.43 


[0.39 0.47] 


0.57 


[0.54 0.61] 


1.0 


[0.99 1.0] 


0.97 


[0.94 1.0] 


a y (Rad s -1 ) 


0.50 


[0.40 0.60] 


0.42 


[0.39 0.46] 


1.0 


[0.99 1.0] 


0.98 


[0.95 1.0] 


a z (Rad s — 1 ) 


0.47 


[0.44 0.49] 


0.39 


[0.38 0.41]] 


1.0 


[0.99 1.0] 


1.01 


[0.99 1.0] 




0.43 


[0.29 0.57] 


0.37 


[0.30 0.43] 


0.13 


[0.13 0.14] 


0.13 


[0.12 0.14] 




0.48 


[0.15 0.81] 


0.37 


[0.28 0.45] 


0.13 


[0.13 0.14] 


0.13 


[0.12 0.14] 


T Lz (s) 


0.39 


[0.33 0.45] 


0.52 


[0.45 0.59] 


0.13 


[0.13 0.14] 


0.13 


[0.12 0.14] 


f n x (Rad 2 s" 1 ) 


0.079 


[0.048 0.11] 


0.12 


[0.093 0.15] 


0.14 


[0.13 0.15] 


0.12 


[0.11 0.13] 


fn B (Rad 2 s^ 1 ) 


0.12 


[0.025 0.21] 


0.066 


[0.048 0.084] 


0.14 


[0.13 0.15] 


0.12 


[0.11 0.13] 


fn, (Rad 2 s" 1 ) 


0.086 


[0.069 0.10] 


0.080 


[0.068 0.092] 


0.14 


[0.13 0.15] 


0.13 


[0.12 0.14] 



Table 5: Calculation of rotational diffusivity, T^. Comparing the laboratory measurements to MC 
simulation show that measurement noise has a slight effect on the calculation of the rotational 
diffusivity. However, the rotational diffusivity is anisotropic, compounded from the anisotropy in 
(Xi and T n% 

5 Discussion 

The hypothesis of this paper is that the rotation of large particles suspended in turbulence satisfies 
an OU process. As mentioned in £J21 the three characteristics of this stochastic process are: expo- 
nentially decaying autocovariance, Gaussian distribution, and stationarity. Other than the noise 
at zero-lag, figures [Th. and [lb show that the exponential decay model is quite effective for particle 
rotation. Based on the discussion of kurtosis above, we conclude that the distribution of angular 
velocities is likely close to Gaussian. We can assume that the rotation process is stationary be- 
cause the driving flow is both stationary and spatially homogeneous. With these three observations 
as support, we conclude that the OU process is an acceptable model for the Lagrangian angular 
velocity of large particles. 

We can use the OU process to determine the time-dependent rotational diffusivity of large 
particles in turbulent flow, T^it). This is given by equation ([2]), and we have obtained af and Th 4 
from the exponential fit to our measurements. Figures Eh. andEb show the time- varying diffusivity 
predicted in this way. When the exponential term diminishes to a negligible value, the Tp^r) 
tends toward a constant value: the Fickian asymptote. The value of the asymptote is calculated 
as = o^Tfi i , and reported in table From these results we see that (for the flow-particle 
combination studied here), particles will diffusively complete one entire rotation cycle in an average 
of roughly 6 minutes (T 27r = ((2vr) 2 /f 0j « 360 s)). 
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Figure 3: Time dependent rotational diffusivity. (a) Spheres and (b) Ellipsoids. Dashed line 
is experimental data (noise removed), solid line is equation ([2]), and dotted line is the Fickian 
asymptote. 



Particle shape appears to have only a minor effect on rotation statistics. The Lagrangian integral 
timescales of angular velocity are statistically identical for ellipsoidal and spherical particles (table 
[3] or [5|) . Due to measurement uncertainty, it is not as conclusive whether th e variance terms are 
statistically similar or different between the two particle shapes. The work of lBellani et al.1 ((2012) 
suggests that the rot ation of particles is dom inated by the large-scales of fluid motion, similar to 



translational motion (jKomasawa et al.l . ll974l ). These results imply that the effect of shape is minor. 



6 Conclusion 

For spherical and ellipsoidal particles suspended in homogeneous, isotropic turbulence, we have 
shown that the particle angular velocity can be modeled as an OU process. This conclusion follows 
from our laboratory measurements showing that the Lagrangian autocovariance of particles' angular 
velocity decays exponentially (figure [T]) . Our measurements also show that the particle rotation 
PDF is close to Gaussian, a conclusion that requires the use of a Monte Carlo simulation to evaluate 
the propagation of measurement noise through our analysis. The results of this demonstrates that 
certain statistics, e.g. the kurtosis, were affected more than others, e.g. the rotational diffusivity. 
Comparing particles of two different shapes, we find that ellipticity does not strongly influence 
rotational diffusion. Both ellipsoids and spheres exhibit similar integral timescales of rotation and 
angular velocity variances. This suggests that the rotational diffusion of large particles in turbulent 
flow can be approximated well without detailed considerations of particle shape. 

We would like to thank Gabriele Bellani and Margaret Byron for support in making the labora- 
tory measurements, Matt Ritter and Audric Collignon for early development of the measurement 
technique, and Jason Lepore, Matt Ford, and Yoram Rubin for insightful discussions. This re- 
search was conducted through undergraduate coursework supported by the department of Civil 
and Environmental Engineering at the University of California, Berkeley. 
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